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ABSTRACT 

New-generation X-ray polarimeters currently under development promise to open a new 
window in the study of high-energy astrophysical sources. Among them, neutron stars ap¬ 
pear particularly suited for polarization measurements. Radiation from the (cooling) surface 
of a neutron star is expected to exhibit a large intrinsic polarization degree due to the star 
strong magnetic field (« G), which influences the plasma opacity in the outermost 

stellar layers. The polarization fraction and polarization angle as measured by an instrument, 
however, do not necessary coincide with the intrinsic ones derived from models of surface 
emission. This is due to the effects of quantum electrodynamics in the highly magnetized 
vacuum around the star (the vacuum polarization) coupled with the rotation of the Stokes pa¬ 
rameters in the plane perpendicular to the line of sight induced by the non-uniform magnetic 
held. Here we revisit the problem and present an efficient method for computing the observed 
polarization fraction and polarization angle in the case of radiation coming from the entire 
surface of a neutron star, accounting for both vacuum polarization and geometrical effects 
due to the extended emitting region. Our approach is fairly general and is illustrated in the 
case of blackbody emission from a neutron star with either a dipolar or a (globally) twisted 
magnetic field. 
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1 INTRODUCTION 


Polarization measurements of radiation coming from astrophysical 
sources helped in improving our knowledge about the physical and 
geometrical properties of a variety of systems, from black holes to 
gamma-ray bursts (e.g. |Trippe|2014| for a review). In this respect, 
neutron stars (NSs) are among the most promising targets for po- 
larimetry due to their strong magnetic field which is expected to 
induce a large degree of polarization of the emitted radiation. 

Radio and optical polarimetry has been already used to derive 
the orientation of the magnetic and rotation axes of radio pulsars 
([Manchester & Taylor||1977[ [Lyne & Manchest^|1988[ see also 
[Pavlov & Zavlin||2000| (. The discovery over the last two decades 
of new classes of X-ray bright, radio-silent NSs with very faint 
(if any) optical counterparts (chiefly the magnetar candidates, e.g. 
[Mereghetti|2008[ [Turolla, Zane & Watts|2015l and the X-ray Dim 


Isolated Neutron Stars, XDINSs, e.g. Turolla 2009[[Kaspi[2010[ ( re¬ 
newed the interest in possible polarization measurements at X-ray 
energies in NS sources. Despite some efforts were made in the past 
to measure polarization in the X-rays, mainly with the OSO-8 and 
INTEGRAL satellites ( [Weisskopf et al.|1978| [Hughes et al.|1984l 
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[Dean et al.[2008[ see also [Kislat et al.[2015[ l, the poor sensitivity of 
past instrumentation did not lead to conclusive results. A new win¬ 
dow opened in the last years, with the advent of new^eneration 
X-ray polarimeters, like XIPE Q IXPE and PRAXyS [j (recently 
selected for the study phase of the ESA M4 and NASA SMEX pro¬ 
grammes respectively), which are based on the photoelectric effect 
and provide a dramatic increase in sensitivity over an energy range 
~ 1-30 keV (see [Bellazzini et al.[2013[ (. X-ray polarimeters derive 
polarization observables by detecting a modulation in the azimuthal 
distribution of events in the focal plane. Actually, while a measure 
of the circular polarization degree is possible in the optical band 
(see e.g. [Wiktorowicz et ar][2015^ , current instruments, based on 
the photoelectric effect or Compton scattering, can only provide 
information about linear polarization ( [Fabiani & Muleri[2014[ l. 

From a theoretical viewpoint, polarization observables (the 

polarization fraction and the polarization angle) are conveniently 
expressed through the Stokes parameters. The comparison between 
the polarization properties of the photons emitted at the source and 
those measured at earth is not straightforward for two main reasons. 
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The first is that the Stokes parameters are defined with respect to a 
given frame, which is in general different for each photon. When 
the Stokes parameters relative to the different photons are added to¬ 
gether, care must be taken to rotate them, so that they are referred to 
the same frame, which coincides with the frame in the focal plane 
of the detector. This effect becomes important every time radiation 
comes from a spatial region endowed with a non-constant magnetic 
field, and will be referred to as “geometrical effect” in the follow¬ 
ing. The second issue, which typically arises in NSs, is related to 
“vacuum polarization”. In the presence of a strong magnetic field, 
quantum electrodynamics (QED) alters the dielectric and magnetic 
properties of the vacuum outside the star, substantially affecting po¬ 
larization ( |Heyl & Shaviv|2002^ . Because of this, (100% linearly 
polarized) photons emitted by the surface will keep their polariza¬ 
tion state up to some distance from the star, as they propagate adia- 
batically. This implies that the degree of polarization and the polar¬ 
ization angle, as measured at infinity, depend also on the extension 
of the “adiabatic region”, which in turn depends on the photon en¬ 
ergy and on the magnetic field. 

The observed polarization properties of radiation from iso¬ 
lated NSs were investigated in the past both in connection with 
the emission from the cooling star surface and the reprocessing of 
photons by magnetospheric electrons through resonant Compton 
scattering, a mechanism which is thought to operate in magnetars. 
IPavlov & Zavlin| ( |2000^ studied the case of thermal emission from 
the entire surface of a NS covered by an atmosphere, without ac¬ 
counting for QED and geometrical effects. A quite complete analy¬ 
sis of the observed polarization properties of surface emission from 
a neutron star has been presented by Heyl, Shaviv & Lloyd|(|200^, 


while |Lai & Ho|p003^ and |van Ade 


sberg & Pema 1 2009 1 focused 


on the role played by the vacuum resonanc^ which occurs in the 
dense atmospheric layers, on the polarization, and may provide a 
direct observational signature of vacuum polarization. The two lat¬ 
ter works were restricted to the case of emission from a small hot 
spot on the NS surface, over which the magnetic field can be treated 
as uniform, therefore no account for rotation of the Stokes param¬ 
eters was required. |Pemandez & Da\^ ( |2011^ and |Tavema et al.] 
( |2014^ have shown that X-ray polarization measurements can pro¬ 
vide independent estimates of the geometrical and physical param¬ 
eters and probe QED effects in the strong field limit in magnetar 


In this paper we re-examine the problem and present a simpli¬ 
fied, efficient method to derive the observed polarization properties 
of radiation emitted from the entire surface of a NS. Our results 
are in agreement with those of |Heyl, Shaviv & Lloyd|(2003| > and 
|van Adelsberg & Pema| ( [2009| l, and our faster approach allows to 

systematically explore the dependence of the polarization observ¬ 
ables on the different geometrical and physical quantities. In partic¬ 
ular, we discuss the difference between the polarization properties 
of the radiation emitted by the star and those measured at earth 
, which is induced by geometrical and QED effects. This aspect, 
which is crucial when one needs to reconstruct the star properties 
from the observed quantities, has not been systematically inves¬ 
tigated in previous works. A complete study based on physically 
consistent models of surface emission is outside the scope of this 
analysis, and we just assume a simple model in which the surface 
emission is a (isotropic) blackbody and the magnetic field is dipo- 


® A Mikheyev-Smimov-Wolfenstein resonance which may induce mode 
conversion in X-ray photons for typical magnetar-like fields {B > 10^^ 

G). 


lar (or a globally twisted dipole field). The outline of the paper is 
as follows. The theoretical framework is introduced in section 
In section [^calculations and results are presented, while section [^ 
contains a discussion about our findings and the conclusions. 


2 THEORETICAL OVERVIEW 

In this section we briefly summarize some basic results about 
the evolution of the polarization state of electromagnetic radiation 
propagating in a strongly magnetized vacuum. Although the con¬ 
siderations we present below are focused on radiation travelling in 
the surroundings of a neutron star, they hold quite in general. 


2.1 Photon polarization in strong magnetic fields 


In the presence of strong magnetic fields photons are linearly polar¬ 
ized in two normal modes: the ordinary mode (O-mode), in which 
the electric field oscillates in the plane of the propagation vec¬ 
tor k and the local magnetic field B, and the extraordinary mode 
(X-mode), in which, instead, the electric field oscillates perpen¬ 
dicularly to both k and B. This holds for photon energies be¬ 
low the electron cyclotron energy {E < E^e = ‘heBjnif.c ~ 


Gnedin & Pavlov 


1974 1 , which implies B > 

10 , 


11.6(5/10^2 G)keV; 

lO^'^ G at X-ray energies, whereas B can be as low as ~ 10^*^ G in 
the optical band. Moreover, the polarization state of photons propa¬ 
gating in vacuo is also influenced by the effects of vacuum polariza¬ 
tion ( |Heyl & Shaviv|2000| [2002] [Harding & Lai|2006[ l. According 
to QED, in fact, photons can temporarily convert into virtual 
pairs. The strong magnetic field polarizes the pairs, modifying the 
dielectric, e, and magnetic permeability, /z, tensors of the vacuum, 
which would coincide with the unit tensor otherwise. 

Fixing a reference frame {x, y, z) with the z-axis along the 
photon propagation direction fc, and the 2 ;-axis perpendicular to 
both k and the local magnetic field B, the evolution of the wave 
electric field is governed by the following system of differential 
equations (see [Fernandez & Davis|2011[ [Taverna et al.|2014[ see 


also |Heyl & Shaviv|2002 for a different, albeit equivalent, formu¬ 
lation) 


dAj; 

dAy 

dz 


ikoS 

ikoS 

~Y~ 


{MA^ + PAy) 
{PA^ + NAy) . 


( 1 ) 


Here A = {Ax, Ay) = {axe~'‘^'° , is the electric 

field complex amplitude, ko = lojc with u) the photon angu¬ 
lar frequency and the adimensional quantities 5, M, N and P 
depend on the (local) magnetic field strength; in particular it is 
5 = {ap/45 tt){B/Bq)^, where ap is the fine structure constant 
and Bq = 4.4 x 10^^ G is the critical magnetic field. As equa¬ 
tions Q show, vacuum polarization induces a change in the electric 
field as the wave propagates; the typical lengthscale over which this 
occurs is = 2/ko6 ~ 100(5/10^'^ G)“2(i<;/1 keV)“^ cm, 
where E — huj. At the same time, the magnetic field changes 
along the photon trajectory, this time over a lengthscale = 
5/|fc • V5| ~ r, where r is the radial distance. Near to the 
star surface it is Ia £b and the direction along which the 
wave electric field oscillates can instantaneously adapts to the vari¬ 
ation of the local magnetic field direction, maintaining the orig¬ 
inal polarization state. In these conditions, the photon is said to 
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propagate adiabatically and in the following we will refer to the 
region in which this occurs as the adiabatic region. However, as 
the photon moves outwards the magnetic field strength decreases 
(B oc + 3 cos^ 9 for a dipole field, where 6 is the mag¬ 

netic colatitude) and £a increases. Since grows more slowly, 
there is an intermediate region in which the wave electric field can 
not promptly follow the variation of the magnetic field any more. 
Finally, in the external region, where I a ^ ^b, the electric field di¬ 
rection freezes, and the polarization modes change as the magnetic 
field direction varies along the photon trajectory. 

The evolution of the polarization modes should be calcu¬ 
lated integrating equations Q from the surface up to infinity (or, 
at least, up to a distance sufficiently large to consider the com¬ 
plex amplitude components and Ay as constants). This has 
been the approch followed by |Heyl, Shaviv & Lloyd| j2003[ see 
also [Fernandez & Davis|20TT]|Tavema et al.|2014) . However, this 
method requires quite long computational times, since numerical 
integration must be carried on along each ray and it is not particu¬ 
larly suited for a systematic study of how the polarization observ¬ 
ables depend on the various physical and geometrical parameters. 
Since the latter is the main goal of the present work, we resort to 
a simpler, approximated treatment in which only the adiabatic re¬ 
gion and the external one are included, and they are divided by a 
sharp edge. To this end we introduce the adiabatic radiu|^ra, de¬ 
fined implicitly by the condition £a = ^b . Assuming a dipole field 
and purely radial photon trajectories, it is Ib = r/3 and hence 
(■Aff = ra/3. Recalling the expression for Ia, it follows that 
ra/7?NS 3.9xl0-'‘(F;/lkeV)-^(Bp/10“ G)-"(7?Ns/ra)“® 
and finally 




IkeV 


1/5 


Rns ■ 


( 2 ) 


where is the stellar radius, Bb is the polar strength of the 
dipole and cos 0 ~ 1 was assumed. The adiabatic radius depends 
on both the photon energy and the star magnetic field: it is larger 
for stars with stronger magnetic field, and, at fixed Bp , it becomes 
smaller for less energetic photons, as shown in Figure[T] 


2.2 Polarized radiative transfer 

A convenient way to describe the polarization properties of the ra¬ 
diation emitted by a source is through the Stokes parameters. With 
reference to the frame (x, y, z) introduced in Section [2T] they are 
related to the complex components of the wave electric field by 

T = AxAj. -\- Ay Ay = -t- 

Q = A^Al - AyA*y =al-al 

(3) 

U = A^A*y + AyA*^ = 2a^ay cos{ipx — >Py) 

V = i {AxAl — AyA*) = 2axay sin((pj: — ipy). 

In the previous equations a star denotes the complex conjugate, X 
is the total intensity of the wave associated to the photon, Q and U 
describe the linear polarization and V the circular polarization. The 
four Stokes parameters satisfy the general relation X^ ^ Q^+U^ + 
V^, the equality holding for 100% polarized radiation. With our 


This same quantity is called the polarization-limiting radius, rpi, in pre¬ 
vious literature (see |Heyl & Shaviv|2002| . 


l0£,n(r^/RM5) 
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Figure 1. Contour plot showing the adiabatic radius (in units of the stel¬ 
lar radius /?ns) a function of the polar magnetic field strength and the 
photon energy. The typical Sp-ranges for different classes of neutron stars, 
the magnetai's, the isolated, thermally emitting NSs (XDINSs), the radio¬ 
pulsars (PSRs) and the central compact objects (CCOs), are also shown. 


current choice of the reference frame, it is = 0 for an X-mode 
photon and ax = 0 for an 0-mode photon. Normalizing the Stokes 
parameters defined above to the intensity X, we can associate to an 
extraordinary/ordinary photon the vectors 


/ Q 
u 


X 





(4) 


where a bar denotes the normalized Stokes parameters. The evolu¬ 
tion of the Stokes parameters mirrors that of the complex compo¬ 
nents of the electric field given in equations 0 (see e.g. [Tavernaj 
jet al.|201^ . Actually, in our hypothesis of 100% linearly polarized 
thermal radiation, the Stokes parameter V is always zero inside the 
adiabatic region. This implies that a circular polarization degree 
can arise only as a consequence of the polarization mode evolu¬ 
tion in the transition between the adiabatic and the external region. 
However, since we do not integrate equations 0 in our model we 
will not discuss V further. We verified that, even accounting for 
the Stokes parameter evolution, as obtained solving equations 0 
across the entire region, the resulting circular polarization fraction 
is very small at optical energies and reaches at most a few percent 
in the X-ray band. 

In order to measure the polarization properties of a given 
source, a polarimeter will collect a large number of photons, each 
characterized by its own set of Stokes parameters. The convenience 
of using the Stokes parameters lies precisely in the fact that they are 
additive (e.g. jRybicki & Lightman|2004| l: the Stokes parameters as¬ 
sociated to the whole collected radiation (i.e. the superposition of 
all the received photons) are equal to the sum of the Stokes pa¬ 
rameters of the single photons. However, care must be taken since 
the quantities in equations ([^ are defined with respect to a precise 
reference frame, {x, y, z), that depends on the direction of the local 
magnetic field B (see & As the magnetic field is in general non- 
uniform across the emission region, its direction at a given point 
will depend on the source magnetic topology. Since the direction of 
the electric field of each photon varies very quickly inside the adia¬ 
batic region, but it is frozen outside, what actually matters is not the 
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Figure 2. Graphical visualization of the different reference frames intro¬ 
duced in the text, k is the common direction of propagation (LOS); (u, 
v) are the fixed, mutually orthogonal axes of the polarimeter; (a:', y') and 
{x", y”) are mutually orthogonal axes of two reference frames relative to 
photons coming from points characterized by the directions B' and B" of 
the magnetic field. The angles a', a" are also indicated. 


B-field direction at the original emission point, but that at the point 
where the photon crosses the adiabatic boundary ra, as pointed out 
by Lai & Ho j2003^ . Let us call {xi,yi,Zi) the reference frame 
in which the Stokes parameters associated to the generic photon 
are defined at the adiabatic radius. While the «i axes are all along 
the same direction (that conicides with the observer line-of-sight, 
LOS), the Xi and yi axes will point, in general, in different direc¬ 
tions for each photon (see Figure]^. To sum correctly the Stokes 
parameters it is necessary to refer them to the same, fixed frame, 
say {u, V, w). This frame can be chosen in such a way to coincide 
with that of the polarimeter, with u and v in the detector plane and 
w along the LOS. 

Each reference frame (xi, j/i, Zi) is rotated with respect to the 
fixed, {u, V, to), frame by an angle a; around the common Zi = w 
axis, where cos Oi = u ■ x;. Under a rotation of the reference frame 
by an angle a;, the Stokes parameters transform as 


7 ; =Ii 


where Nx (No ) is the number of extraordinary (ordinary) photons, 
N = Nx + No, and we used equations 


2.3 Polarization observables 


The polarization state of the detected radiation can be described in 
terms of two observable^ the linear polarization fraction Hl and 
the polarization angle Xp defined as 


Hl 


^/wnp 

I 


1 

Xp = - arctan 



(7) 


The linear polarization fraction is not, in general, equivalent to the 


ratio I Nx—No \ /N (as previously noticed by Heyl, Shaviv & Lloyd 


|2003| l. This would happen only if all the angles Oi were the same, 
i.e. when the magnetic field is uniform across the emitting region 
(as in the case considered by|Lai & Ho|2003|and[van Adelsberg &| 
|Pema|2009| of radiation coming from a small hot spot on the NS 
surface). In fact, denoting with ap the common value and substitut¬ 
ing expressions in the first of equations l|^, one obtains 

y/{Nx — Wo)^ cos2(2ao) + {No — Wx)^ sin^(2ao) 


Hl = 


N 


\Nx-No\ 

N 


( 8 ) 


Under the same hypothesis the polarization angle, given by the sec¬ 
ond of equations is directly related to the angle qq 


1 

Xp = 2 arctan 


(Nq - Nx) sin(2Qo) 
{Nx — No) cos(2ao) 


—ao . 


(9) 


Hence, the polarization fraction gives direct information about the 
intrinsic degree of polarization of the radiation (i.e. that at the 
source) only for a constant rotation angle ap. Under the same con¬ 
ditions, the polarization angle provides the direction of the (uni¬ 
form) magnetic field of the source in the plane of the sky. 

On the contrary, if the iJ-field is non-uniform (e.g. for emis¬ 
sion coming from the entire surface of a NS endowed with a dipole 
field), a will vary according to the magnetic field direction at the 
point where the photon crosses the adiabatic radius. In this case 
equations 0 and 0 give 


Hl = 


N 


N 


+ 2 cos(2ai — 2Q;k) 


-f2y^ cos(2aj — 2Q!r) 

- 2cos(2ai - 2aj) 


1/2 


( 10 ) 


Qi = QiCOs(2ai) +Wisin(2Q!i) (5) 

Ui = Ui cos(2Q;i) — Qi sin(2Q!i). 

Since photons emitted by the star surface, or, more generally, inside 
the adiabatic region, are 100% polarized either in the X or O mode 
(i.e. Qi = ±1 and Ui — 0), the Stokes parameters of the radiation 
collected at infinity are 

Q = Qi = cos(2ai) - cos(2aj) 

( 6 ) 

u = Ui = E ^°1 sm(2«j) - ES sin(2ai) 


where i, k = 1,..., Nx and j, r = 1,..., No, while the polarization 
angle results in 


Xp = 


1 

- arctan 
2 


E.^"isin(2«i)-E.gisin(2aj) - 
E^Ei cos(2ai) - cos(2aj). 


( 11 ) 


So, in the general case both Hl and Xp depend on the distribution 
of the angles Oi, which, in turn, is determined by the geometry of 
the magnetic field. 


® As mentioned earlier, circular polarization is not considered in the 
present work. 
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3 POLARIZATION OF SURFACE EMISSION EROM 
NEUTRON STARS 

In this section we present quantitative results for the polarization 
observables in the case of surface (thermal) emission from a neu¬ 
tron star endowed with an axially-symmetric magnetic field, either 
a dipole or a (globally) twisted dipole, the latter often used to de¬ 
scribe the magnetosphere of magnetars (see|Thompson, Lyutikov| 
|& Kulkarni|2002l l. 


3.1 The a-distribution 

Let us introduce a reference frame (X, Y, Z) with the Z-axis in the 
direction of the LOS (unit vector t), X in the plane of £ and the 
star spin axis (unit vector f2) and Y — x X. The geometry is 
shown in Figure]^, where x is the angle between the spin axis and 
the LOS, and ^ is the angle between the spin axis and the magnetic 
dipole axis (unit vector bdtp)- It is f2 = (— sin x, 0, cos x) while, 
having introduced the polar angles rj and ^ that fix the direction 
of bdip with respect to t (see again Figure |^), one has bdip = 
(sin rj cos sin rj sin cos rf). The angles rj and (( are related to x 
and ^ by 


magnetic field in this frame are 

(17) 

where 6 is the magnetic colatitude. Then, the cartesian components 
B = (Bp, Bq,Bt) can be calculated making use of expressions 
l |A3^ in Appendix [A| However, since all the calculations to derive 
the analytical expression of the angle a (see equation |16[) are in 
the LOS reference frame, the X, Y and Z components of B are 
needed. They can be obtained through a change of basis as 

Bx = Bppx + BqQx -I- Bt{bdip)x 

By = BpPY + BqPY + Bt(bdip)v (18) 

Bz = BpPz + BqPz + Bt{biip)z , 

where the components of p, and q in the {X, Y, Z) frame are given 
in Appendix]^ while those of bdip are given by equation ^13[ l. 

Substituting the expressions l |18| l in the first of equations 0 , 
the components of the unit vector x in the LOS reference frame are 


fBj, + 


-By \ 

Bx . (19) 

0 / 


,pol _ 





cos r) = cos X cos ^ -I- sin x sin ^ cos 7 

c (12) 

^ cos E — cos X cos 77 
cos C = " , 

sin X sin 77 

where 7 is the rotational phase. Using the previous expressions, the 
components of the unit vector bdip in the frame [X, Y, Z) become 

( sin X cos 5 — cos x sin ^ cos 7 
sin ^ sin 7 

cos X cos ^ -I- sin x sin ^ cos 7 

According to the discussion in §2.2| the axes u and v of the po- 
larimeter frame can be chosen as any pair of orthogonal directions 
in the XY plane. In general it is 




/ cos tp \ 

1 

^ — sin p) \ 


u = X^ 

sin 7 /) 

, v=Y={ 

cosip 

. (14) 


\ 0 

\ 

V 0 / 



where 7 /) is the angle that the u axis makes with the X axis. The 
axes X and y of the reference frame ( 2 :, y, z), that change for each 
photon, are defined once the magnetic field geometry is fixed as 


_ £x B 

^ \l X B\ 


y — X X . (15) 


The angle a by which the photon frame {x, y, z) has to be rotated 
to coincide with the polarimeter frame (m, v, w) is then simply ob¬ 
tained taking the scalar product of u with x 


cos a = u ■ X . 


(16) 


The indetermination in the sign of a is resolved looking at the sign 
of 77 ■ a;; if the latter is positive the rotation is by an angle —a (i.e. 
sin a = —x /1 — cos^ a). 

Since we need to consider photons only from the boundary of 
the adiabatic region outwards, B in the first of equations (HI is the 
stellar magnetic field calculated at (see equation]^. Actually, it 
is more convenient to express the magnetic field components in a 
reference frame (p, g, t), with the t axis along bdip and p, q two 
mutually orthogonal directions in the plane perpendicular to t (see 
Figure]^). For a dipole, in particular, the polar components of the 


Bx and By clearly depend on the angles x. C srid Ihe phase 7 
through the unit vectors of the (p, q, t) reference frame, given by 
equations l [B2l l, l |B^ and Moreover, they depend also on the 
magnetic colatitude and azimuth (9 and (p that fix the point where 
the magnetic field is calculated on the adiabatic surface) through 
the components Bp, Bq and Bt given by equations jA3[ > and ^17[ l. 
Actually, the angles 6 and (j) depend in turn on x, ^ and 7. To make 
this dependence explicit, let us consider Figure that shows, in 
the LOS reference frame, the path of a photon emitted from a point 
of the surface characterized by the polar angles 0 s and <I?s, up to 
the point where it crosses the adiabatic boundary, characterized by 
the angles 0 and 3>. Observing the star at infinity, one collects only 
photons that travel along vectors k = (0,0, k) parallel to the LOS 
£. The modulus of each vector k is fixed by the condition: 

ro -I- fc = ra , (20) 

where ro = i?Ns(sin 0s cos T>s, sin 0s sin T>s, cos 0s) is the po¬ 
sition vector of the surface point from which the photon has been 
emitted and ra = ra(sin 0 cos sin 0 sin cos 0 ) is the po¬ 
sition vector of the point where the photon crosses the adiabatic 
boundary. Taking the norm of both the sides of equation l |20[ l and 
solving for k, the only acceptable solution is 

k= \Jrl- 77^3 sin^ 0 s - Rns cos 0s ( 21 ) 

and, substituting this result again in equation po^ , one obtains: 

( i?NS sin 0s cos 'Fs \ 

77ns sin 0 s sin -Fs j , ( 22 ) 

\/ri- Bf s sin^ 0 s / 

where the distance ra of the adiabatic boundary is given by equa¬ 
tion iD- From simple geometrical considerations (see again Figure 
[^), it follows that 

cos 9 = bdip • — , (23) 

T’a 

while the cosine of the angle (j> can be obtained as 

cos = p ■ ra , (24) 

where r^^ is the unit vector of the projection of ra orthogonal to 
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Figure 3. The two reference frames used in the calculation of the a angle. Left: the (X, y, Z) frame with the Z axis in the direction of the LOS the X axis 
in the plane where H is the star spin axis, and the Y axis perpendicular to both X and Z. Right: the (p, q, t) reference frame with the t axis along the 
star magnetic axis 6dip, the p and q two mutually orthogonal axes in the plane perpendicular to bdip (see Appendix[B]for more details). The angles C’ 

6 and (f) are also shown. 



Figure 4. The path of a photon emitted by a point on the star surface with 
polar angles ©s and that crosses the adiabatic boundary in a point of 
polar coordinates © and 


6 dip. The complete expressions of cos 6 and cos 0 are given in Ap¬ 
pendix]^ 

Finally, substituting into equation gives the distribution 

of a 


cos a = 


Bx sin Ip 


By cos Ip 


(25) 


which is a function of the angles x. the phase 7, the photon 
energy E and Bp (through the adiabatic radius r^, see equation 
1 ^, the polar angles 0s and $s that fix the point on the surface 
from which the photons were emitted and the angle ip by which 
the polarimeter frame is rotated wrt the LOS one. In the following 
we take ip = 0, i.e. the u (v) axis coincides with the X (Y) axis, 
although the generalization to other values is straightforward. 


3.2 Numerical implementation 

In order to calculate the polarization fraction Hl and the polariza¬ 
tion angle Xp, we use the ray-tracing code developed by|Zane &| 
|Turolla| ( |2006^ , with the addition of a specific module for the eval¬ 
uation of the a angle distribution and of the Stokes parameters. 
QED effects are included as described in section]^ The code takes 
also into account the effects due to the strong gravity on photon 
propagation (relativistic ray-bending) and on the stellar magnetic 
field. For a dipole field (see equations|17|l, the latter are given by 


pGR _ f p 

-^r — Jdip-^r 


Bf^ = gdipBg 

B°^ = B^ = 0, 
where 

^x{x + 2) 

— ) 
l-xj ’ 


/dip — 


ln(l — a;) -I- 


— ,/i _ ^ ( _ 9 r,. 


(26) 


(27) 


with X = Rs/r’, Rb = 2 GMns jf? is the Schwarzschild radius and 
Mns is the stellar mass (see |Page & Sarmiento|1996l >. 

The expressions for the total Stokes parameters Q and U given 
in equations can be easily generalized to a continuous photon 
distribution by replacing the sums with integrals over the visible 


VB'i + 
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Figure 5. Contour plots in the energy-phase plane of the polarization fraction (top row) and the polarization angle (bottom row) for a neutron star with 
BP = 10^^ G, i?NS = 10 km and mass Mns = 1-4 Mq. The inclination of the magnetic axis with respect to the spin axis is fixed to ^ = 5°, while the 
angle between the spin axis and the LOS is x = 0° (left column), 15° (middle column) and 30° (right column). Seed photons are 100% polarized in the 
X-mode. 


part of the star surface 


pOtir n 1 

Fq= d$s / du^(nx — no) cos(2a) 

Jo Jo (28) 

r 27 r rl 

Fu= d$s / du^(no — nx) sin(2a), 

Jo Jo 

where nx (no) is the photon intensity in the extraordinary (ordi¬ 
nary) mode and Fq and Fu are the “fluxes” of the Stokes parame¬ 

ters (see |Pavlov & Zavlin|2000| l. In general, nx and no depend on 
the photon energy E and direction, and on the position on the star 
surface of the emission point. The integration variable u = sin 0 
is related to 0s by the integral (see |Zane & Turolla|2006| and ref¬ 
erences therein): 


0 = 



dn sin 0s 

[(1 — a:) /4 — (1 — 2vx) sin^ 0s] 


(29) 


that accounts for ray-bending and reduces to 0 = 0s in the limit 
a; —>■ 0 (when the effects of general relativity can be neglected). 
The total photon flux is obtained in a similar way 

^27r n 1 

Fi= d^s / dw^(nx-I-no) . (30) 

Jo Jo 

For the sake of simplicity, we assume in the following that 
radiation is emitted by the cooling star surface with an isotropic 
blackbody distribution. The photon intensity is then 


nx,o 


2 E'^ 

e^/kT _ ^ 


(31) 


where T is the local surface temperature. In order to model the 


surface thermal distribution and to avoid a vanishing temperature 
at the equator, here we adopt a variant of the standard tempera¬ 
ture distribution for a core-centred dipole field (e.g. |Page||1995l >, 
Titt) = max(Tp| cos Te), where i9 is the angle between the 
local normal and B, Tp and Te are the temperature at the pole and 
at the equator, respectively; in the following we take Tp = 150 eV 
and Te = 100 eV. The polarization degree of the radiation emitted 
at the surface is fixed specifying the ratio po = tix/(nx + no), 
|nx - no|/(nx -f no) = |2po - 1|- 


3.3 Results 

The polarization observables IIl and Xp can be computed recall¬ 
ing the definitions given in equation 0 and using the expressions 
we just derived for the Stokes parameters, equations ( |28^ and j30| l, 
together with the distribution of cos a given in equation l|2^. All 
results presented in this section refer to a neutron star with mass 
ATns = 1 .4M0 and radius i?NS = 10 km. Thermal photons 
are assumed to be 100% polarized in one of the two modes, i.e. 
Po = 0,1; essentially we consider all the photons as extraordi¬ 
nary, unless explicitly stated otherwise (see e.g. the discussion in 
[Fernandez & Davis|2011[ [Tavema et al.|2014l l. This is the choice 
which produces the most unfavourable conditions to detect the de¬ 
polarizing effects of vacuum polarization and geometry on the po¬ 
larization observables. 

Figure]^ shows the polarization fraction and the polarization 
angle as functions of the photon energy and the rotational phase for 
different values of the inclination x of the LOS wrt the star spin 
axis. The magnetic axis is at an angle ^ = 5° with respect to the 
spin axis (i.e. the NS is a nearly aligned rotator) and Bp = 10^® 
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Figure 6. Polarization angle as a function of the rotational phase at a fixed 
energy (E = 0.02 keV), for x = 90° and different values of g: 0° (red), 
30° (orange), 45° (green), 60° (light blue), 75° (blue) and 90° (violet). 
The solid (dashed) lines correspond to seed photons 100% polarized in the 
X-mode (0-mode). The values of i?NS. ^NS ^nd Bp are the same as in 
Figure]^ 


G. The effects produced by the frame rotation (induced by the non- 
uniform B-field) are quite dramatic, as it is evident from the polar¬ 
ization fraction (top row). In particular, for x = 0 (top left panel) 
IIl is almost everywhere far from unity, the value expected from 
the intrinsic degree of polarization, |nx — no|/(nx + no) = 1, 
and it becomes ~ 0.9 only at i? ~ 10 keV. By increasing the LOS 
inclination (x = 15°, top middle panel), the polarization fraction 
reaches unity for photon energies > 1 keV, while at lower energies 
it is substantially smaller (between ~ 0.1 and ~ 0.8). Only when x 
becomes sufficiently large (x = 30°, top right panel) IIl is unity, 
except at low energies (~ 1-10 eV), where the polarization fraction 
drops to about 0.6 in some phase intervals. 

The bottom row of Figurej^shows Xp for the same three simu¬ 
lations. Contrary of what happens for the polarization fraction, the 
polarization angle does not depend on the energy and exhibits an 
oscillatory behavior as a function of the rotational phase around a 
value of 90°. The amplitude of the oscillations depends on the geo¬ 
metrical angles and, for some combinations of x and Xp sweeps 
the entire range [0°, 180°] through a discontinuity, or “jump”. This 
is clearly seen in the bottom left panel of Figure]^ where x = 0 , 
while Xp is in between ~ 70°-110° (bottom middle panel) and 
~ 80°-100° (bottom right panel) for x = 15° and x = 30°, re¬ 
spectively. This is further illustrated in Figure]^ which shows the 
polarization angle as a function of the rotational phase at a single 
energy {E = 0.02 keV), x = 90° and different values of ^ for 
radiation 100% polarized in the X-mode (solid lines) and in the O- 
mode (dashed lines). The amplitude of the oscillation vanishes in 
the case of an aligned rotator seen equator-on (x = 90°, ^ = 0°) 
and increases for increasing ^ until the “jump” appears for ^ = 90° 
The average value of Xp, instead, does not change with x and ^ 
and is fixed by the polarization mode of the seed photons: it is 90° 
for X-mode photons and 0° for O-mode ones[^ It should be noted, 

® The curves for ^ = 90° in Figure|^are box-like; the sloping lines are an 
artifact introduced by the finite resolution of the phase grid. 

^ Actually the mean value is the same even if photons are not all polarized 
in the same mode; since the Stokes parameters for 0- and X-mode photons 


however, that the mean value of the polarization angle is not univo- 
cally associated to the two photon modes, since it depends on the 
choice of the fixed reference frame {u,v,w), i.e. on the angle t/) 
introduced in j j3.1| If, for instance, t/) = 90° (so that the u axis co¬ 
incides with the Y axis of the LOS reference frame), the situation 
depicted in Figure|^is reversed, with the polarization angle for X- 
mode photons oscillating around 0° and that for the O-mode ones 
around 90°. Of course, different choices of the ip angle do not af¬ 
fect the polarization degree IIl, the amplitude of the oscillations of 
Xp and the shift of 90° between the mean values of xp for X-mode 
and O-mode photons. 

The behaviour of the phase-averaged polarization fraction as 
a function of the angles x> ^ is shown in Figure |7] for two values 
of the energy, S = 2 eV (optical) and E = 0.3 keV (X-rays). The 
right panel illustrates the variation of the semi-amplitude of xp . As 
already noted by |Femandez & Davis| ( |20lT| , the amplitude is 180° 
for 5 < X when the phase-averaged polarization degree attains its 
minimum value (see 0- 

The effects of varying the magnetic field strength are illus¬ 
trated in Figure]^ where x = 15°, C = 5° and Bp = 10^^ G (left 
panel). Bp — lO'^® G (middle panel; this is the same case shown 
in Figure]^ and Bp = 10^'* G (right panel). Again, changes are 
mostly in the polarization fraction IIl (top row). Overall, the polar¬ 
ization fraction is smaller when the magnetic field is lower (top left 
panel), and increases for increasing Bp, reaching values ~ 1 (i.e. 
the intrinsic polarization degree) in almost the entire energy range 
for Bp — 10^'* G (top right panel), see 0 On the contrary, the 
polarization angle (bottom row) does not change much, exhibiting 
an oscillation between ~ 70° and ~ 110° at all the values of Bp. 

Finally, Figure|^illustrates the effects on the polarization ob¬ 
servables induced by the presence of a toroidal field component. 
The right column shows the phase-energy contour plot of the po¬ 
larization fraction (topjsanel) and a phase plot of the polarization 
angle at a fixed energ) j(bottom panel) for a globally twisted dipole 
field. The left column shows for comparison the same quantities for 
a pure dipole with the same Bp = 10^® G. The twisted magnetic 
field was evaluated using the analytical approximation by lPavan"^ 
|aL| ( |2009] see expressions in their Appendix A), with a twist angle 
A())n-s — 50°. Since relativistic corrections are unavailable for 
a twisted field, they were not applied also to the dipole we show 
for comparison, whereas ray bending is still considered in both the 
cases. The effects of the twist on the polarization fraction are quite 
modest and the variation of IIl with photon energy and rotational 
phase is nearly the same as in the pure dipole case. The only differ¬ 
ence is in a slight overall decrease in the polarization degree. 

The twist of the external field affects much more the polariza¬ 
tion angle, as it can be seen from the bottom row of Figure]^ The 
net effect is an overall asymmetry of the oscillations: Xp sweeps a 
larger angle in a half-period with respect to the purely dipolar case; 
this effect increases with the twist angle Acpt^s, as already noticed 
by |Femandez & DavisH201 f) . 


have opposite signs and the polarization observables are obtained by sum¬ 
ming the Stokes parameters over all photons, the mean value of Xp reflects 
the polarization mode which dominates. 

® The twisted field actually introduces a dependence of Xp on the photon 
energy. For the values of the twist angle we consider, however, this depen¬ 
dence is quite small. 
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Figure 7. Contour plots for the phase-averaged polarization fraction at optical (2 eV, left panel) and X-ray (0.3 keV, middle panel) energies, and of the semi¬ 
amplitude of the oscillations of the polarization angle (right panel), as functions of x ^nd The values of ^ns -^P l^e same as in Figure 
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Figure 8. Same as in Figure|^for x = 15°. ^ = 5° and three different values of the magnetic field: fjp = 10^^ (left column), 10^® (middle column) and 
10^"^ G (right colmun). 


4 DISCUSSION AND CONCLUSIONS 


In this paper we reconsidered the problem of the relation between 
the intrinsic and observed polarization properties in the case of 
surface emission from a neutron star. Our work extends previous 
investigations l Heyl, Shaviv & Lloyd|2003l |Lai & Ho|2003[ |van| 


[Adelsberg & Perna|2009^ by providing the polarization observables 
for a large set of physical and geometrical parameters (i.e. the an¬ 
gles X and the photon energy and the magnetic field strength 
and topology). Our treatment includes both “geometrical” effects, 
due to the rotation of the Stokes parameters which is needed when 
the magnetic field is not constant across the emitting region, and 
“vacuum polarization” ( |Heyl & Shaviv|2002| ); in order to make a 
full exploration of the parameter space possible, an approximated 


treatment of QED was used. This resulted in much shorter compu¬ 
tational times, without loosing significant physical accuracy, as the 
comparison with available results shows ( |Heyl, Shaviv & Lloyd| 
|2003| >. We checked that a typical run required about 100 minutes 
integrating equations Q and only few tens of seconds using our 
approximation. Moreover, our approach allows to better disentan¬ 
gle the effects of QED and those due to the rotation of the Stokes 
parameters on the polarization signals. We stress that our main goal 
was not to compute polarization observables for a precise, physical 
model of surface emission, but to systematically illustrate the role 
of these two effects in making polarization patterns different from 
those of the original radiation. To this end we assumed a quite sim¬ 
ple picture in which the star magnetic field is a dipole and each 
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Figure 9. Polarization observables for the cases of a pure dipolar magnetic field (left column) and a globally twisted dipole with twist angle A(?)n_s — 50° 
(right column). Top row: polarization fraction in the energy-phase plane for X = 15° and ^ = 5°. Bottom row: polarization angle as a function of the 
rotational phase for a fixed photon energy (E = 0.02 keV), x = 90° and g = 0° (red), 30° (orange), 45° (green), 60° (light blue), 75° (blue) and 90° 
(violet). All the plots are obtained for seed photons 100% polarized in the X-mode and values of ^NS and Bp as in Figure]^ 


surface patch emits a (isotropic) blackbody spectrum at the local 
temperature. Our results can be easily generalized to other mag¬ 
netic configurations (the case of a globally twisted field is actually 
discussed here) and to different surface emission models. A (com¬ 
parative) analysis of the polarization observables for emission from 
an atmosphere (van Adelsberg & Lai 2006 and references therein) 
or from a condensed surface (Potekhin et al.|2012[ and references 
therein) will be the subject of a future paper (Gonzalez Caniulef et 
al., in preparation). 


A crucial point in assessing the measured polarization prop¬ 
erties is that the polarization state of a photon propagating in a 
magnetized vacuum (either ordinary or extraordinary) is strictly re¬ 
lated to the choice of the reference frame. In fact, the photon po¬ 
larization mode is defined only with respect to the plane fixed by 
the wavevector and the local magnetic field. This means that the 
Stokes parameters of each photon are in general referred to differ¬ 
ent frames, with the two axes orthogonal to the direction of propa¬ 
gation tailored on the direction of the local B-field ( §2.2[ ). However, 
photons are collected in the focal plane of an instrument, where a 
reference direction has been a priori introduced. This means that 
in order to obtain the polarization observables relative to the pho¬ 
tons received by the instrument in a given exposure time, the Stokes 
parameters of each photon must be transformed to the polarimeter 
reference frame, through a rotation in the plane orthogonal to the 
line-of-sight by an angle a, which depends on the magnetic field 


and viewing geometry, on the photon energy and on the position of 
the point from which photons were emitted ( §3.1^ . 

The effects induced by rotation of the reference frames com¬ 
pound those of vacuum polarization. According to QED, in fact, 
photons maintain their initial polarization state within the adiabatic 
region, while the polarization freezes at a larger distance (see & 
Despite the transition between the adiabatic and the outer zones is 
smooth, we assumed that there is a sharp boundary at the adiabatic 
radius r^, (see equation]^. This enabled us to treat the photons as 
if they were emitted at r^., as far as their polarization state is con¬ 
cerned. This implies that the distribution of the a angles, by which 
each frame has to be rotated, is actually determined by the magnetic 
field at the adiabatic radius, and hence depends also on the distance 
Ta from the star surface. 

Simulations of polarization measurements ( §3.3^ clearly show 
that, because of the combined effects of frame rotation and QED, 
the measured polarization fraction can be very different from the 
intrinsic value, i.e. that of the radiation emitted at the surface. The 
differences appear to depend firstly on the viewing geometry, i.e. 
on the angles x and ^ which give the inclination of the line-of- 
sight and of the dipole axis with respect to the star spin axis. As 
shown in Figure]^ the polarization dramatically decreases at all 
rotational phases for x — C — 0. For x not too close to Hl 
has a minimum at the phase 7 where the magnetic axis bdip lies 
in the plane of the rotation axis and the LOS (either 7 = 0 or 27r 
with our choice of the reference frame). This behaviour confirms 
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the results of |Heyl, Shaviv & Lloyd| ( [2003^ and it is entirely due 
to the non-constant magnetic field across the emitting region: when 
a region of the star near to the magnetic poles is into view, the 
projection of B in the plane orthogonal to the LOS is essentially 
radial, so that a can take values in the entire range [0, 27r]. Since the 
Stokes parameters for the whole radiation are obtained integrating 
the rotated Stokes parameters of single photons over the part of 
the star in view (see equations |28|l, the polarization degree has a 
minimum when the angle between bdip and the LOS is minimum 
(and equal to \x — ^|, see the first of equations |12[>. On the other 
hand, the fact that IIl does not change with rotational phase in the 
case shown in the top-left panel of Figurej^ is precisely due to the 
fact that this is a nearly aligned rotator viewed along the rotational 
axis. 


In agreement with the results by |Heyl, Shaviv & Lloyd|p003| l, 
we found that the behavior of IIl is also sensitive to the loca¬ 
tion of the adiabatic radius. From the top rows of Figures]^ and 
it can be seen that the linear polarization fraction increases with 
the photon energy and the polar strength of the B-field: this re¬ 
flects the dependence of on In fact, as equation (25 

shows, cos a depends on the magnetic co-latitude and azimuth, Q 
and (j>, through Bx , By ; the two latter angles contain the factor 
Bxs/fa (see equations Cl and C3 '. So, in the limit ra S> Rns (at 


least for axisymmetric magnetic field topology), a remains nearly 
constant as the emission point changes on the star surface, imply¬ 
ing that the polarization fraction can be indeed approximated with 
|nx — no\/(nx^ -f no), as equation ^ shows. Heyl, Shaviv & 


|Lloyd|p003| l explained this behavior as due to the fact that, in this 
limit, QED birefringence aligns the photon polarization angles. Ac¬ 
tually, the weaker depolarization when ra S> Rns is due chiefly to 
the rotation of the Stokes parameters, vacuum polarization entering 
only implicitly through the dependence of the angle q on ra. Be¬ 
cause of the dependence of ra on E and Bp, this approximation 
becomes better the larger the photon energy and the stronger the 
polar B-field. Instead, the closer to the star surface the adiabatic 
limit, the smaller the overall measured polarization degree, the lat¬ 
ter becoming vanishingly small if no adiabatic region is accounted 
for. So, the main conclusion is the more point-like the star is seen 
by an observer at the adiabatic boundary, the closer the measured 
IIl is to the intrinsic linear polarization degree. A similar effect 
was noted by |Heyl, Shaviv & Lloyd|f2003) in connection with the 
variation of the stellar radius. 

On the other hand, the polarization angle exhibits quite a dif¬ 
ferent behaviour. As Figures and show, Xp do^s not change 
significantly with ra, since it does not depend on E and Bp. This 
is because the factors Rns /fa within cos a tend to cancel out tak¬ 
ing the ratio U/Q which defines Xp (equation]^. The fact that Xp 
keeps oscillating even when the measured polarization fraction is 
much smaller than the intrinsic one (see bottom-left panels of Fig- 
ures[^and[^ is a consequence of the frame rotation and not of QED 
effects. The polarization angle depends quite strongly, instead, on 
the geometrical angles x and ^ (see e.g. the right panel of Eigure 
1^. The polarization swing generally increases for decreasing x at 
fixed as shown in Eigure In particular, Xp sweeps the entire 
range [0,180°] when the region close to the magnetic pole is al¬ 
ways in view during the star rotation (bottom left panel), while the 
swing gets smaller for values of x and ^ such that the polar region 
enters into view only at certain rotational phases. On the other hand, 
the oscillation amplitude in general grows for increasing ^ at fixed 
X- This behavior appears to be related again to the a-angle distri¬ 
bution, and provides an explanation for the correlation between the 
swing by 180° of the polarization angle and the low phase-averaged 


polarization fraction at x < already noticed by |Ee r nandez &| 
|Davis|j201 1| see also |Wagner & Seifert|2000) . In fact, the regions 
where the polarization angle spans the widest range correspond to 
those in which at least one among the Stokes parameters Q and U 
takes all the values between —1 and 1. Consequently, the averaged 
polarization fraction, obtained by summing the Stokes parameters 
over a rotational cycle, turns out to be very small, as shown in Eig¬ 
ure [T] 

Phase-resolved polarization angle measurements, together 
with the information given by the linear polarization fraction, can 
help in understanding which polarization mode is the dominant one 
in the detected radiation. We showed in Figure that the mean 
value of Xp depends on the mode in which the majority of pho¬ 
tons are polarized. However, it is also related to the orientation of 
the {u,v) axes in the polarimeter plane (the t/; angle, see ]3.1 1 , 
which are fixed by the instrument design. In particular, the mean 
values of Xp for X- and O-mode photons are always displaced by 
90°, but they are 90° and 0° respectively (as in the case in Figure 
1^ only if tp = 0. Hence, a measurement of the polarization angle 
alone fails in telling which is the prevailing polarization mode. The 
problem can be solved if also a phase-resolved measurement of the 
linear polarization fraction is available. In this case, since Hl has 
a minimum when bdip intercepts the ft — £ plane (see above), it 
could indeed be possible to individuate the direction of the X axis 
on the plane of the sky. This allows to derive the angle ip and to 
remove the inherent ambiguity in the measurement of Xp- 

Polarization observables can also provide information on the 
source geometry, i.e. the inclination of the LOS and of the mag¬ 
netic axis wrt the rotation axis. In fact, as discussed earlier on, both 
the polarization fraction and the polarization angle strongly depend 
on the angles x and As already shown in |Tavema et al.| ( [2014| ), 
if phase-resolved polarization signals are available, a simultaneous 
fit of Hl and Xp (possibly supplemented by that of the flux) allows 
to unequivocally derive the values of x and Q On the contrary, this 
is in general not possible starting from phase-averaged measure¬ 
ments. The phase-averaged polarization fraction is largely degen¬ 
erate with respect to the two angles, as clearly shown in Figure]^ 
and, since the phase average polarization angle is constant in large 
regions of the x ~ C plane, its measure is of no avail in pinpointing 
X and Q 

The effects of a different magnetic field topology on the 
polarization observables were assessed in the illustrative case of 
globally-twisted dipolar magnetic field[^ The presence of a toroidal 
component in the external magnetic field slightly changes the be¬ 
haviour of linear polarization fraction (see Figure]^. In a twisted 
field, depolarization induced by the frame rotation is a bit stronger. 
This is due to the fact that Btwist > Bdip at any given position, 
because the toroidal component is roughly of the same order of the 
poloidal one, while the r-dependence is about the same for the two 
magnetic configurations. As a consequence the adiabatic boundary 
moves a bit closer to the surface if Bp and the photon energy are the 
same. A twisted field influences the polarization angle most, pro¬ 
ducing a strong asymmetry in the swing and a weak dependence 
on the energy (mainly at optical energies), as already discussed by 
[Fernandez & Davi^P011[ > and jTavema et al.| ( [2014t . The fact that 
the polarization angle is more sensitive to QED effects for a twisted 


® We focused here only on surface emission, the interactions of photons 
with magnetospheric currents, chiefly through resonant cyclotron scatter¬ 
ing, were ignored. 



























12 R. Taverna et al. 


magnetosphere than for a purely dipolar field, provides a strong sig¬ 
nature of vacuum polarization effects (see |Taverna et al.|2014^ . 

Our analysis further demonstrates the need to properly ac¬ 
count for QED and frame rotation effects in evaluating the observed 
polarization properties of radiation emitted by a neutron star. This 
is of particular relevance in relation to recently proposed X-ray po- 
larimetry missions, which will certainly select neutron star sources 
as primary targets. 
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APPENDIX A: CARTESIAN COMPONENTS OE B 

The cartesian compononents of the magnetic field B in the ref¬ 
erence frame {p, q, t) can be obtained from its polar components 
= {Br, Be, B^) given in equation l |l7| using the following 
expression 

B = (Bp, B^,Bt) = (BP°' • • qP°', 5^°* • 4^°*), 

(Al) 

where and are the unit vectors relative to the (p, q, t) 

frame expressed in polar components 

— p ■ {f, 6, (j)) = (sin 6 cos (j), cos 9 cos 0, — sin (jj) 

= q ■ = (sin 9 sin 0, cos 9 sin cos (j>) (A2) 

^pol _ ^ ^ ^ 

and the angles 9 and 4> are the magnetic colatitude and azimuth, 
respectively (see Figure]^). Upon substituting expressions l |A2| > 
into equation ([ATJ, one finally obtains 

Bp = sin 9 cos 4>Br + cos 9 cos (j)Be — sin (j)B,i, 

Bq = sin 9 sin (f>Br + cos 9 sin (f>B 0 + cos cjiB^ (A3) 

Bt = cos 9Br — sin 9Be . 


APPENDIX B: MAGNETIC REEERENCE ERAME 


The projection of bdip, given by equation ( |13f , orthogonal to the 
spin axis f2 (see ] 3.1 1 , in the LOS reference frame {X, Y, Z) is 


6dip — (bdip • 

|bdip — (bdip • $T)$T| 


— cos X cos 7 
sin 7 

sin X cos 7 


(Bl) 


m is an unit vector corotating with the star around the spin axis. 
The projection of m perpendicular to bdip fixes the p axis of the 
bdip reference frame (p, q, t). Its expression in the {X, Y, Z) frame 
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is given by 


P = 


m- [m- bdip)bdip 
m — (m- bdip)bdip| 

( — sin X sin ^ — cos x cos ^ cos 7 
cos ^ sin 7 

sin X cos ^ cos 7 — cos x sin ^ 


(B2) 


Finally, the unit vector defining the q axis, in the (X, Y, Z) refer¬ 
ence frame, is given by the vector product between bdip and p 


Q — bdip X p — 


- COS X sin 7 
— COS 7 
sin X sin 7 


(B3) 


APPENDIX C: COMPLETE EXPRESSIONS EOR cos 6 

AND COS(j) 

Substituting the components of bdip (equation|131 and of ra (equa- 
tion| 22 ^ into equation P3^ one obtains 

cos 8 — sin 0s ( cos >l>s sin x cos ^ -|- sin $3 sin ^ sin 7 
ra V 

— cos d>s cos X sin ^ cos 7 


1 f Rns . 

1 — -sm 0 > 


V 


sin 0s ) ( cos X cos ^ -I- sin x sin ^ cos 7 

(Cl) 


To calculate the complete expression of cos (j) we use equation j24| l, 
where the components of the unit vector p are given in equation 
I and 


(C2) 


sin $s cos ^ sin 7 — cos $s sin x sin ^ 
— cos <l>s cos X cos ^ cos 7 
ri - {Rns sin0s)2 


^_L _ fa (ra ' bdip)bdip 

|ra - (ra • bdiplbdipT 
Putting all together, one obtains 
-Rns sin 0 s 


cos (/) = 


+ 


Ta sin 6 


rj sin^ 9 


sin X cos 5 cos 7 — cos x sin ^ j . 

(C3) 








